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I.  INTRODUCTION 


In  the  treatment  of  the  response  of  aircraft  to  atmospheric 
turbulence,  a common  assumption  has  been  to  consider  that  the  air- 
craft system  is  represented  by  a rigid  body  with  the  single  degree 
of  freedom  of  vertical  motion  only.  Many  design  studies  and  de- 
sign charts  have  been  based  on  this  assumption;  for  example, 
references  1,  2 and  3.  The  consideration  of  several  degrees  of 
freedom,  including  vertical  motion,  pitch,  and  flexible  modes,  has 
generally  been  limited  to  isolated  treatments  of  specific  aircraft, 
such  as  shown  in  reference  4.  It  should  be  noted  that  the  results 
for  the  response  parameters  A and  N0  obtained  from  a multi- 
degree-of- freedom  analysis  may  differ  considerably  from  the  results 
of  a single  degree -of- freedom  analysis;  the  results  should  there- 
fore not  be  used  interchangeably,  as  is  often  done. 

The  general  intent  of  multimode  treatments  is  to  bring  out 
changes  in  response  or  amplification  effects  due  to  flexibility. 

Such  response  analyses  thus  tacitly  concentrate  on  higher  frequency 
effects.  It  Is  found,  however,  that  the  low  frequency  effects  may 
be  even  more  significant  than  the  high  frequency  effects.  Consider 
figure  1.  Flexibility  may  cause  a rather  large  change  in  the  air- 
plane transfer  function  at  the  higher  frequencies  (as  compared  with 
the  rigid  body  result),  but  this  change  is  weighted  by  rather  small 
values  of  the  Input  spectrum.  By  contrast,  the  change  in  the 
transfer  function  at  low  frequencies  due  to  including  pitch  may 
appear  small,  but  this  change  is  weighted  by  very  large  values  of 
the  input  spectrum.  The  resulting  effect  which  shows  up  In  the 
output  spectrum  of  Including  or  not  including  pitch  may  therefore 
be  greater  than  flexibility  effects. 

The  effect  of  including  pitch  may  be  seen  somewhat  by  physical 
reasoning.  If  vertical  motion  only  is  considered,  the  airplane  tends 
to  follow  or  go  up  and  down  with  the  gusts  In  the  low  frequency 
range.  With  pitch  included,  however,  the  airplane  acts  as  a 
weathercock  at  low  frequencies,  with  the  result  that  there  is 
virtually  no  vertical  motion  response.  In  effect,  the  inclusion  of 
the  degree-o f- freedom  of  pitch  with  vertical  motion  causes  the 
response  spectrum  to  shift  to  higher  frequencies. 

The  basic  point  being  made  is  that  not  enough  attention  has 
been  focused  on  the  low  frequency  end  of  the  gust  response  spectrum. 
Reference  5 gives  a treatment  involving  vertical  motion  and  pitch. 

The  analysis  is  based,  however,  on  a stability  derivative  approach, 
and  it  is  not  known  how  well  the  analysis  applies  at  intermediate 
frequencies,  where  most  of  the  gust  response  power  seems  to  appear. 

The  purpose  of  this  proposed  Investigation  is  thus  to  make  a 
systematic  study  of  gust  response,  based  on  spectral  techniques 
wherein  both  the  degrees-of- freedom  of  vertical  motion  and  pitch 
are  included.  The  aim  is  to  derive  parametric-type  charts  which 
allow  for  the  easy  evaluation  of  the  response  parameters  A and  N . 


pj.W-rPPB  i|  ii.iip,. 


A.  Aerodynamic  Loading 

The  aerodynamic  forces  used  herein  are  found  in  accordance 
with  the  analysis  given  in  reference  6.  The  development  in  this 
reference  is  in  terms  of  the  average  downwash  over  a finite  span 
interval  due  to  a concentrated  or  point  aerodynamic  load.  An 
equivalent  procedure  is  to  establish  the  downwash  at  a point  due 
to  a uniform  line  load  of  finite  length.  This  equivalent  procedure 
is  used  herein. 


The  aerodynamic  loading  considered  in  this  study  is  shown  in 
figure  2;  this  figure  also  serves  to  indicate  some  of  the  notation 
used.  It  is  to  be  noted  that  an  aim  in  this  study  was  to  keep  the 
number  of  aerodynamic  loads  to  a minimum,  yet  still  of  sufficient 
generality  so  as  to  represent  overall  loads  and  moments  on  the 
airplane  in  a realistic  manner  (the  aim  is  to  establish  the  air- 
plane dynamics  reasonably  well,  not  to  obtain  precise  airloads  over 
the  wing).  The  wing  loading  is  seen  to  be  represented  by  two  line 
loads,  while  a single  line  load  is  used  for  the  tail.  Later 
examples  will  show  that  this  rather  "crude''  representation  of  the 
aerodynamics  gives  fairly  accurate  loads.  The  following  consider- 
ation serves  to  indicate  that  two  line  loads  on  the  wing  should  be 
adequate  for  taking  into  account  the  gust  frequency  components  of 
concern.  Consider  a gust  frequency  component  with  wavelength  A ; 
spatial  frequency  then  follows  as 


n - 

These  relations  lead  in  turn  to 
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where  k ■ . Past  gust  studies  (reference  2)  show  that  the 

largest  value  of  k of  concern  is  around  k ■ .5  ; this  value 
leads  to  the  shortest  wavelength  of  concern  of 
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Relative  to  this  shortest  wavelength,  two  control  or  downwash  points 
in  the  chordwise  direction  would  appear  as  in  the  following  sketch. 


Effectively,  one  wavelength  of  the  shortest  gust  component  is 
represented  by  12  points,  which  is  a fairly  good  approximation. 
All  lower  k values,  or  longer  wavelengths,  where  there  is  more 
response  power,  are  of  course  represented  even  more  accurately. 


The  control  or  downwash  points,  two  on  the  wing,  one  on  the 
tail,  that  are  used  in  the  subsequent  development  are  also  shown 
in  figure  2.  Definition  of  two  aspect-ratio-type  terms  are 
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If  the  wing  is  assumed  to  behave  in  the  fashion  of  an  elliptical 
wing,  some  additional  parameters  of  convenient  use  may  be  es- 
tablished, such  as  the  following: 
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Use  will  be  made  of  these  relations  in  later  examples. 
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B.  Aerodynamic  Forces 

The  aerodynamic  forces  that  are  of  concern  in  this  study  are 
derived  in  this  section.  The  strength  of  the  wake  potential  due 
to  a uniform  line  load  of  intensity  p , see  figure  3(a),  is  given 
by 

_ lux 

+ <!> 

Since  this  line  load  extends  from  y«-jtoy»^,  the  total 
load  is  given  by 

P « pX  (2) 

For  present  application  purposes,  we  need  the  downwash  values 
along  the  x-axis,  both  ahead  of  and  behind  the  line  load.  To  es- 
tablish these  downwash  values,  use  is  made  of  the  potential  for  a 
unit  dipole  source,  given  by 


Mir  (x  +y*+z  )J/* 

The  downwash  in  the  z « 0 plane  that  is  associated  with  this 
dipole  source,  for  regions  outside  the  origin,  is  given  by 


(x*+y*) 


By  equations  (1)  and  (3),  the  downwash  ahead  of  the  line  load 
(for  y « 0)  is  given  by  (see  figure  3(b)) 
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This  equation  may  be  converted  to  the  following  result 


where 
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in  which 


For  the  downwash  behind  (downstream)  the  line  load  we  first 
establish  the  potential  in  the  field  due  to  the  strength  given  by 
equation  (1).  By  equations  (1)  and  (3)  the  field  potential  for 
y ■ 0 is  given  by  (see  figure  3(c)) 


k/i  * 

f f 


[(x-5)2  ♦ y2  ♦ x2]' 


It  is  convenient  to  change  the  variable  of  integration  in  the  x 
direction,  and  then  rewrite  the  equation  so  that  it  assumes  the 
form 


1(0  Xr  *>eo  X/2  00  CO  oo-i 
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where 


f - - e 

IT 


i^jX 


(x^+y^+z*) 


By  transforms  918.5  and  867  in  reference  7,  the  first  term  within 
the  brackets  of  equation  (8)  can  be  integrated,  the  result  being 


simply  2e 


; equation  (8)  thus  appears 
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The  downwash  in  the  z « 0 plane  indicated  by  this  equation  is 
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where 


(x*  + yc) 


Equation  (10)  may  be  integrated  further*  and  then  reduced  to  the 
following  result 
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where 
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cos  ks 


(i 2 ) ds 


(lib) 


and  where  k and  s are  as  before. 

Equations  (6)  and  (11)  are  the  basic  equations  used  herein 
for  the  evaluation  of  the  downwash  that  results  from  the  aerodynamic 
loads  on  the  wing  and  the  tail.  For  later  application,  it  is  con- 
venient to  Introduce  a simplified  but  descriptive  notation  of  these 
equations,  both  in  application  to  the  wing  and  the  tail.  For  the 
wing,  let 


wm  “ fpUcI^mn  + i®mn^ 


where  the  notation  is  interpreted  as  the  downwash  at  point  m due 

to  a load  at  point  n . The  definitions  for  a and  8 are 

mn  mn 

as  follows:  For  a downwash  point  ahead  of  the  line  load  we  have. 
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from  equation  (6) 


a * jg  •.  e (C,  - iSn ) 
mn  pnm  ' 1 X7 


while  for  a downwash  point  behind  the  line  load, equation  (li)glvaa 
the  definition 

-w,  ■ -'1!:S;co ♦ C1  • 1S1>  'w 

It  ic  to  be  noted  that  the  values  of  s used  in  equation  (13)  is 
based  on  the  absolute  value  of  the  distance  the  downwash  point  is 
ahead  of  the  line  load;  likewise  s in  equation  (14)  is  based  on 
the  absolute  value  of  the  distance  the  downwash  point  is  behind 
the  line  load. 

In  the  application  of  eouations  (6)  and  (11)  to  the  tail, 
consideration  must  oe  given  to  the  fact  that  > and  c mav  be 
different.  For  the  tail 

2x 

t ct 
uc. 
kt  * w~ 

and  to  express  the  k values  in  terms  of  a common  value,  we  write 

kt r-  k (15) 

In  this  study  downwash  points  on  the  tail  ahead  of  the  line  load 
are  not  considered.  For  downwash  points  behind  the  tail  line 
load,  equations  (11)  and  (15)  give 
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Steady  State  Cheok.~  As  a check  on  the  accuracy  that  might  be 
expected  from  the  assumed  line  load  distribution  used  herein, 
application  was  made  to  a wing  in  a steady  state  constant  angle  of 
attack  condition.  Two  cases  were  studied:  one  with  a single  line 
load;  one  with  two  line  loads,  as  shown  in  the  following  central 
cross  sections  of  the  wing. 
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CASE  I 


1* 


W1  *P 2 w2 


CASE  II 


3c/8-| 

- 5c/8  - 
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Results  obtained  by  applying  equations  (6)  and  (11)  are  as  follows. 
For  Case  I,  P was  found  to  be  given  by 


*pU2Sctr 


1 + /I  + a4 

where  a denotes  angle  of  attack. 
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Consider  a wing  with  aspect  ratio  A = 8 , then 


a a IF  * 8 * **-9348 

P in  turn  evaluates  to 

P - . 8177ttpU2Scx0 

For  an  elliptical  wing  of  A « 8 , and  with  the  use  of  the  factor 
a 

for  aspect  ratio  correction  to  the  slope  of  the  lift  curve, 

P would  evaluate  to 

P - .8TrpU2SaQ 

The  result  obtained  from  using  only  a single  finite  length  line 
load  is  seen  to  be  only  about  25?  different  than  the  corresponding 
exact  solution  for  an  elliptic  wing. 

For  Case  II,  the  following  results  were  found: 

P^  ■ . 6l6il'TrpU2Sa0 
P2  » .2013TTpU2Sao 
lPn  « P ■ .8l77TrpU2Sa0 

The  result  is  seen  to  be  the  same  as  for  Case  I, 

C.  Equations  of  Motion 

In  reference  to  figure  2,  the  equations  for  vertical  and 


pitching  motion  of  the  airplane  are 

mz  « ?1  + P2  + p3  (17) 

mr2$  - (e  + §)P1  - (|£  - e)P2  - etP3  (18) 

The  vertical  displacements  of  the  control  points  w.  , w?  , and 
w^  are  given  by 

z1  - z + (e  - !>)♦  (19) 

z2  ■ z + (e  - |^)0  (20) 

c. 

z^  ■ z - (e^  + (21) 
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The  downwash  at  the  control  points  Is  due  to  three  sources:  the 
vertical  velocity  movement  as  can  be  found  from  equations  (19), 
(20),  and  (21),  the  instantaneous  angle  of  attack,  and  the 
sinusoidal  gust  which  is  assumed  to  flow  over  the  airplane  (for 
purposes  of  determining  the  frequency  response  function  due  to 
sinusoidal  gust  encounter).  Consideration  of  these  three  sources 
yields  the  downwash  values 


w^  « 4 + (e  - -jj)4>  - U$  - wQ  (22) 

^OJC 

w.,  * z + (e  - |£)$  - U<f>  - wQe  2U  (23) 

a to 

• ^4-  • “"iff 

* z - (e^  + )4  - U<£  “ wQe  (2^) 

where  wQ  is  the  magnitude  of  the  sinusoidal  gust,  and 

el=et+e"§+r‘ 


We  note  that  these  equations  take  into  account  the  lag  in  the  time 
that  the  sinusoidal  gust  is  felt  at  points  2 and  3 relative  to 
point  1. 

By  equations  (12),  (13),  ( 14 ) and  (16)  the  downwash  at  the 
control  points,  when  expressed  in  terms  of  the  loads  P,  , P2  , 
and  , may  also  be  written  as 

W1  * yC(cxll  + l0ll)Pl  + (a12  + i012)P23  (25) 

w2  * y a21  + ^21^P1  + ^ On  + ^ll^P2-)  (26) 

w3  * yC(«3i  + 1^3i)pi  + (a32  + l832)P2  + tjT(a33  + le33)P33 

(27) 

where  y « TrpUS  ; note  c^2  “ a11  , and  B22  * • We  note  that 

the  downwash  on  the  tail  due  to  the  wing  loads  is  taken  into  account; 
the  downwash  on  the  wing  due  to  the  tail  load  is  ignored,  however, 
on  the  basis  that  this  effect  should  be  small. 

We  now  set  equations  (22),  (23),  and  (2H)  equal  to  equations 
(25),  (26),  and  (27).  The  three  equations  thus  formed,  together 
with  equations  (17)  and  (18)  then  form  a set  of  5 equations  with  5 
unknowns,  z , $ » pi  » P2  and  P3  • The  equations  so  found  are 

given  in  Table  I (in  a nondimensional  representation,  and  with  wQ 
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set  equal  to  unity  to  apply  to  a unit  sinusoidal  gust).  The 
application  of  equations  (12),  (13),  ( 14 ) and  (16)  to  the  three 
specific  control  points  shown  In  figure  2 Involved  the  choice  of 
s values  as  follows : 

For  s 


all  & eil 


1 

2 


a12  & 612 


1 

2 


a21  & S21 


a _ , & R 

31  31 


a32  & 032 


a33  & 633 


1 + 2(£  + + ft 

¥ c c c 


-I 


+ 2(-  + —•)  + — £■ 
vc  c c 


It  should  be  noted  that  the  equations  given  in  Table  I repre- 
sent a rather  unique  formulation  of  the  gust  response  problem.  We 
do  not  establish  the  loads  that  are  due  to  airplane  motion  and  due 
to  gust  encounter,  and  then  feed  these  loads  into  the  dynamical 
equations  of  motion  to  establish  t frequency  response,  as  is 
usually  the  conventional  approach.  Here,  we  let  the  loads  be  un- 
known quantities  just  as  the  variables  z and  <p  . All  aero- 
dynamic effects  - nonsteady  lift  effects,  spanwlse  induction 
effects,  tail  downwash,  and  the  lag  of  the  gusts  in  traveling  from 
the  wing  to  the  tail  - are  all  automatically  taken  into  account. 
For  solution,  all  we  need  to  solve^for  is  uz  , from  which  the 

vertical  load  factor  follows  as  . Direct  determination  of 

S 


the  loads 


P1  , » and 


is  completely  bypassed.  We  note, 


however,  that  the  equations  allow  for  the  ready  determination  of 
the  aerodynamic  loads  on  the  airplane,  if  desired.  All  we  have  to 
do  is  consider  the  last  3 equations;  solution  of  these  equations 
for  P1  , P2  , and  P^  gives  the  loads  that  are  due  to  vertical 

motion  ujz  , or  due  to  pitch  $ , or  due  to  gust  encounter  w . 


We  might  note  how  the  extension  to  two  degrees  of  freedom  has 
greatly  Increased  the  number  of  airplane  parameters  involved.  For 
the  case  of  the  single  degree  of  freedom  of  vertical  motion,  the 
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airplane  can  be  describe;!  by  one  parameter  only,  the  mass  parameter 
P (reference  2).  For  the  case  of  the  two  degrees  of  freedom  of 
vertical  motion  and  pitch,  at  least  seven  parameters  are  required 
to  describe  the  airplane,  namely: 


where  r is  the  radius  of  gyration  for  pitching  inertia.  This 
large  number  of  parameters  essentially  precludes  the  generation 
of  generalized  parametric  charts. 
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K 
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D.  Evaluation  of  the  Response  Spectrum 

The  evaluation  of  the  response  spectrum  follows  readily  from 
the  equations  in  Table  I as  follows.  First,  consider  the  determi- 
nation of  the  frequency  response  function  for  c.g.  vertical 
acceleration  due  to  sinusoidal  gust  encounter.  Through  use  of  the 
equations  given  in  Table  I,  the  resulting  c.g.  acceleration  is 
given  as 


g 

The  output  spectrum  for  An  is  then 

*An  " !Anl?*w  * 1^1  *w  (28) 


is  the  gust  input  spectrum.  The  rms  value  for  An 

\ r f 0 2 2 "1 1/2 

v -LJ  iv  vy  (29) 

0 

where  kc  is  some  cut-off  frequency  (to  be  discussed  sub- 
sequently). We  wish  to  rearrange  equation  (29)  so  that  it  has  a 
form  similar  to  that  used  in  past  gust  studies;  specifically,  the 
following  form  is  sought 


where  <p 

w 

follows  as 


a 


An 


ttqSU 

W 


Ko 


1 


where 


U_  K 
eg  y al 


TTpegS 


(30) 


and  where  K is  in  the  nature  of  an  alleviation  factor,  and  a, 
is  some  rms  gu3t  input  value.  When  the  rearrangement  is  made,  * 


I 

'I 
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we  find  that 


[p.*-r 


(3D 


where 

f1  * 4y2k2  [ o>z  ( 2 

The  value  of  k which  is  associated  with  the  number  of  up- 
ward crossings  per  second  according  to  the  equation 


N = k 
o ire  o 


follows  as 


r /Vf.  % cue  1 1/2 

k ■ i o; 

° ?e  ♦ 

L / fi  -f  dk  J 

0 ai 

For  use  in  computational  studies,  the  value  of  kc  used  will  be 
the  rule-of-thumb  value  that  is  discussed  in  reference  8,  namely 

i tt 
kc  3 A 

where  A as  used  in  this  definition  refers  to  wing  aspect  ratio. 

The  gust  input  spectrum  that  is  adopted  herein  is  given  by 
the  equation 


Vk) 


a2  1 ♦ f(l.339  It)* 

w 


1 + (l. 


2L  \2 
339  f*  kj 


This  equation  may  be  rearranged  so  as  to  yield  the  result 


♦w<“> 


(i^[lt  1(1.339 

1 + (l.  339  k)2J11/6 


0 


H 


where 


(34) 


A plot  of  equation  (33)  is  shown  in  figure  4.  We  note  that  the 
rearrangement  was  made  so  that  all  the  spectra  pass  through  the 
same  values  at  high  k , regardless  of  the  value  of  2L/c  . This 
alteration  in  spectrum  definition  is  made  because  of  the  belief 
that  results  obtained  may  be  fairly  insensitive  to  cr  independent 
of  the  2L/c  value  used. 
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III.  RESULTS 


i 


A.  Check  with  One-Degree-of-Freedom  Case 


As  a partial  check  of  the  equations  given  in  Table  I,  solution 
was  made  for  the  case  where  the  pitching  motion  was  set  equal  to 
zero  (set  <}>  - 0 and  strike  out  the  second  equation);  thus,  a 
solution  for  vertical  motion  response  alone  is  obtained.  Typical 
results  for  the  f^  function  are  shown  in  figure  5.  The  open 
points  are  the  results  for  y * 20  and  y » 100  as  given  by  the 
single-degree-of- freedom  treatment  of  reference  2.  It  is  noted 
that  good  agreement  is  obtained  throughout.  The  slight  difference 
in  results  is  probably  due  to  the  fact  that  the  analysis  presented 
herein  takes  into  account  lag  in  lift  and  gust  penetration  effects 
in  different  but  perhaps  better  fashion.  The  main  point  of  the 
figure  is  that  it  serves  to  indicate  that  the  basic  equations  are 
correct. 


B.  Results  for  Basic  Cases  Studied 

To  establish  an  understanding  of  how  pitch  influences  the  air- 
plane response  due  to  gusts,  four  cases  covering  a range  in  the 
basic  parameters  were  studied.  The  following  table  indicates  the 
parameters  used  in  these  four  basic  cases. 


Case 

A 

a 

at 

e 

c 

_t 

c 

ol  o 

At 

X" 

S 

37 

r 

c 

I 

6 

3.701 

2.056 

0 

3.3 

.6 

1/3 

5 

1 

II 

6 

3.701 

2.056 

.15 

3.15 

.6 

1/3 

5 

1 

III 

10 

6.169 

3.427 

0 

3.3 

.6 

1/3 

5 

1 

IV 

10 

6.169 

3.427 

.15 

3.15 

.6 

1/3 

5 

1 

e 


These  cases  were  chosen  mainly  to  bring  out  the  effect  of  aspect 
ratio  A , and  of  the  importance  of  the  c.g.  position  as  given  by 
e 

c ' 


One  of  the  most  significant  results  of  the  present  study  is 
shown  in  figure  6.  In  this  figure  the  f^  function,  which  repre- 
sents the  square  of  the  amplitude  of  the  frequency  response 
function,  is  given  for  the  airplane  with  two  degrees  of  freedom 
(vertical  motion  and  pitch)  and  for  the  airplane  with  vertical 
motion  only.  Two  significant  points  are  seen.  First,  there  is  a 
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very  pronounced  peak  for  the  two-degree-of-freedom  airplane, 
associated  with  the  short  period  mode.  Second,  the  frequency 
response  function  for  the  vertical  motion  only  case  Is  very  large 
at  low  frequencies,  whereas  the  function  is  negligible  at  low  fre- 
quencies for  the  two-degree-of-freedom  case,  due  to  the  weather- 
cocking  of  the  airplane,  as  mentioned  in  the  introduction.  Gust- 
power  at  the  low  frequencies  is  thus  significantly  lower  for  the 

two-degree-of-freedom  airplane  as  compared  with  the  airplane  with  \ 

the  single  degree  of  vertical  motion  only. 

Figure  7 compares  how  the  response  power  for  c.g.  vertical 
acceleration  is  distributed  with  frequency  for  the  vertical  motion 
only  case  and  for  the  airplane  with  pitch  and  vertical  motion.  In 
this  figure  k is  P-1-otte^  against  log  k ; this  type  plot  is 
significant  because  it  shows  how  the  response  power,  as  given  by 

the  area  under  the  curves  as  seen,  is  distributed  with  frequency.  ; 

Two  points  of  significance  are  noted.  First,  the  major  response 
power  for  the  two-degree-of-freedom  case  appears  at  much  higher 

frequencies  in  comparison  to  the  vertical  motion  only  case.  Second,  ‘ 

the  response  for  the  vertical  motion  only  case  is  strongly  dependent  { 

on  the  ratio;  for  the  two-degree-of-freedom  case,  the  response  j 

c 2L  I 

appears  to  be  rather  insensitive  to  the  value  of  ~ . 


r 


Values  of  K and  kQ  obtained  for  Case  III  are  shown  as  a 

function  of  the  mass  parameter  p In  figures  8 and  9.  In  the  study 
of  these  figures  we  should  keep  in  mind  that  practical  values  of 
2L 

— appear  to  be  100  or  greater  (suppose  L = 1000  ft,  and  c ■ 10  ft ; 
^ 2 L pT 

then  — * 200).  Thus,  if  we  discount  the  — * 50  curve,  and  con- 
c c 

slder  u value  less  than  about  100,  which  represents  the  range  of 

practical  interest,  we  see  that  the  results  for  K and  k0  can  be 

pT 

represented  essentially  by  single  curves  independent  of  ~ . This 

observation  Is  in  marked  contrast  to  results  that  are  found  for  the 
vertical  motion  case  (reference  2).  The  fact  that  the  consideration 
of  both  pitch  and  vertical  motion  seems  to  give  gust  response  re- 
sults which  are  essentially  independent  of  the  turbulent  scale  L 
Is  a significant  finding  of  this  study.  The  results  suggest  that 
we  can  eliminate  L , which  has  been  the  subject  of  much  controversy 
over  the  past  several  years,  as  one  of  the  parameters  to  be  con- 
sidered In  the  design  for  gusts. 


Results  found  for  Cases  I,  II  and  IV  are  very  similar  trend- 
wise  to  the  results  shown  in  figures  8 and  9,  and  are  thus  not 
presented.  Examination  of  the  results  for  all  four  cases  Indicates 


that  the  K or  k curve  for 


2L 

- uux’ve  i ur  - — 
O C 

purposes  as  the  single  curve  which 
2L 

values  of  — greater  than  100  . 


■ 200  can  be  taken  for  practical 
represents  the  results  for  all 
On  this  basis,  figures  10  and  11 


t f***1  —"t 


. . TT*"  **il 


16 


PT 

give  the  £=■  = 200  results  for  K and  Ko  for  all  the 

basic  cases  studied.  These  figures  indicate  that  gust  response  is 
dependent  to  some  extent  on  wing  aspect  ratio  and  on  c.g.  position, 
as  might  be  expected. 


results  for 


for  all  the 


C.  Results  for  Parameter  Variations 

The  basic  cases  studied  in  the  previous  section  covered  mainly 
aspect  ratio  and  c.g.  position  effects.  In  this  section  results 
are  given  for  a systematic  variation  of  other  airplane  configuration 
parameters.  The  main  intent  is  to  gain  an  idea  of  the  sensitivity 
of  gust  response  to  the  various  parameters  (we  see  that  the  mass 
parameter  y is  still  a very  significant  parameter).  Six  cases, 
as  listed  in  the  following  table,  were  studied. 


Case 


Nominal  3.3  1 

A 3.63  1 

B 2.97  1 

C 3.3  1 

D 3. 3 1 

E 3.3  .9 

F 3.3  1.1 

Common  parameters  to  all  six  cases  were 


10  * I 


Cases 


A and  B represent  a 10  percent  increase  and  a 10  percent  de- 
crease in  the  tail  moment  arm  relative  to  the  nominal  case.  In 
cases  C and  D the  ratio  of  wing  area  to  tail  area  is  de- 
creased and  increased  10  percent.  Cases  E and  P represent  a 
10  percent  decrease  and  a 10  percent  increase  in  the  radius  of 
gyration  in  pitch. 


Results  for 


are  shown  for  the  six  cases  in 


figures  12  and  13,  as  obtained  using  ~ « 200  . There  is  some 

c 

variation  in  the  results,  but,  in  general,  the  variation  is  small. 
In  the  range  of  y of  practical  interest,  say  up  to  y of  about 
60  , it  would  appear  feasible  to  represent  the  results  by  a single 
curve,  with  errors  less  than  5 percent. 
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D.  Effects  of  Downwash  on  the  Tail 


To  study  how  important  tail  downwash  effects  are,  two  cases 
were  studied,  one  with  tail  downwash  effects  included,  one  with 
the  effects  removed.  In  reference  to  the  equations  given  in 
Table  I,  the  ct^i  , Sgi  , 01^2  * and  e32  terms  are  associated  with 

tail  downwash  effects.  To  neglect  the  influence  of  downwash  on  the 
tail,  we  simply  set  these  terms  equal  to  zero.  Results  for  kf^$w 

as  obtained  for  basic  Case  III  are  shown  in  figure  14.  It  is  seen 
that  the  inclusion  of  downwash  causes  the  peak  to  shift  to  a lower 
frequency,  and,  as  Judged  by  the  width  of  the  peak,  causes  an  in- 
crease in  damping  in  pitch.  Results  for  K are  shown  in  figure  15. 
The  neglect  of  downwash  is  seen  to  increase  the  gust  response. 


E.  Comparison  with  Vertical  Motion  Only  Case 

It  is  of  interest  to  see  how  the  results  of  the  present  study 
compare  with  previously  established  results  for  an  airplane  having 
the  single  degree  of  freedom  of  vertical  motion  only.  If  we 
examine  equation  (33),  it  can  be  shown  that  the  K values  given  in 
this  report  are  related  to  the  K values  given  in  reference  2 by 
the  relation  ¥ 


This  realtion  yields 


■Mir)3 


.1531 

.1216 

.0965 

,0766 


By  these  conversion  numbers,  values  of 


were  established  for 


Case  III  from  figure  8.  Results  are  shown  in  figure  16  in  compari- 
son to  the  results  for  the  vertical  motion  only  case  as  taken  from 

2L 

reference  2.  It  is  noted  that  for  — » 50  , the  K . values  for 

C 9 

the  vertical  motion  case  are  greater  than  the  values  for  the 

two-degree-of- freedom  airplane  at  low  y , but  become  significantly 

less  at  high  y values.  For  ■ 400  , the  values  for  the 

vertical  motion  case  are  significantly  greater  over  the  entire 
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U range  shown.  This  figure  shows  that  the  A values  computed  for 
the  two-degree -of- freedom  case  can  be  considerably  different  than 
the  A values  for  the  case  of  vertical  motion  alone.  The  impli- 
cation is  that,  either  in  the  design  for  gust  or  in  the  reduction 
of  gust  data,  consistent  procedures  for  establishing  A and  Nq 

should  be  used.  For  example,  it  is  improper  to  compare  deduced 
gust  data  obtained  by  one  investigator  using  one  procedure  to 
compute  A with  gust  data  obtained  by  another  investigator  who 
used  a different  procedure  for  establishing  A . 

In  figure  17  results  for  kQ  are  compared.  The  solid  curve 
applies  to  the  two-degree-of-freedom  case  and  represents  an 
"average”  of  the  results  given  in  figure  11;  the  dashed  curve 
applying  to  the  vertical  motion  only  case  is  taken  from  reference  8 
(A  = 8).  It  is  seen,  in  contrast  to  the  findings  for  A , that  the 
results  for  k0  for  the  two  different  cases  are  in  fairly  good 
agreement.  This  is  not  surprising,  since  the  response  spectrum  for 
the  two  cases  is  essentially  the  same  at  high  frequencies  (see 
figure  7),  and  since  k depends  on  a k^  weighting  of  the 
response  spectrum. 

Because  the  spectral  function  used  herein  was  chosen  so  as  to 

2L 

be  independent  of  -=■  for  large  values  of  k (see  equation  (33) 

V 

and  figure  k ) , the  values  of  K found  in  this  study  are  different 
by  an  order  of  magnitude  than  the  corresponding  K values  found  in 
other  studies.  It  is  of  interest  to  see  what  level  of  gust  design 
velocities  should  be  used  with  the  K of  the  present  report  to 
yield  acceleration  levels  which  are  similar  to  those  found  by  the 
discrete-gust  design  approach.  For  the  discrete-gust  design 
approach  we  have 


Analogous  to  this  equation,  and  on  the  basis  of  the  form  suggested 
by  equation  (30),  we  write 

An  ■ ^§SV 

where  K refers  to  the  values  found  in  this  report,  and  is 

the  associated  design  gust  velocity.  We  note  that  in  this  report 
the  mass  parameter  p is  defined  in  terms  of  the  theoretical  lift 
curve  slope  of  2tr  j in  this  section  we  designate  this  parameter 
with  a subscript,  or 

Wl  « W 

1 TTpCgS 

to  distinguish  it  from  the  mass  parameter  used  in  the  discrete-gust 
approach,  namely 
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As  an  example,  consider  an  airplane  with  a ■ 5 and  y * 30  . The 

value  y^  would  then  be  y^  ■ y ■ 23.9  • For  y " 30  the 

value  of  Kg  Is  .75  (reference  2),  while  figure  10  yields  a 

value  of  K of  aroun'd  4.5  (using  y ■ 23.9).  If  we  equate  the  two 

jn  equations,  and  use  Ud  - 50  , we  find 

5 x .75  x 50  - 2*  x 4.5  x Ux 

which  yields  U,  » 6.63  . Thus,  aj propriate  gust  design  values  for 

use  with  K values  as  derived  in  this  report  are  in  the  neighbor- 
hood of  6.5  to  7.0  fps. 

It  is  of  interest  also  to  show  how  the  results  of  the  present 
report  are  converted  so  as  to  yield  the  structural  response  para- 
meter Ar  as  used  in  reference  3*  In  reference  3,  o^n  is  given 

a.  ■ A 0 
An  r w 

If  this  equation  and  equation  (30)  are  combined,  the  following 
result  for  A is  found 

A .iLiL.fi 

r v1  ow 

where  y-^  13  defined  as  used  in  this  section.  Equation  (34)  yields 


where  n is  defined  as  used  in  equation  (35).  Thus,  Ap  becomes 

A - 3L  IL  n 

r eg  y2 

For  computational  purposes,  we  select  a specific  value  of  n ; the 

value  suggested  is  n ■ .12  (which  corresponds  to  a — value  near 
100);  the  specific  Ar  equation  is  thus  c 

Ar--12^  <36, 


i 


w 


By  contrast,  the  value  of  Ar  as  given  in  references  2 and  3 
(single  degree  of  freedom  of  vertical  motion)  is  given  by 


A » 
r 


u_5t 

eg  u 


(37) 


An  example  application  is  the  following.  Suppose  that 


v1  * 23.9  , or  u - 30 

C * 20  ft 
U = 500  fps 


The  results  of  this  report  show  that  K is  approximately  4.5  for 
“ 23.9  ; references  2 or  3 Indicate  ^ ■ .0205  for  u * 30 
(with  * ICO).  Equation  (36)  then  yields 

Ar  - .0175 

while  equation  (37)  yields 

Ar  * .0159 

In  this  specific  example,  the  inclusion  of  pitch  with  vertical 
motion  is  seen  to  increase  the  response  parameter  A as  compared 
to  the  result  for  vertical  motion  alone.  r 


IV.  CONCLUDING  REMARKS 


The  inclusion  of  the  pitching  degree  of  freedom  with  vertical 
motion  is  found  to  considerably  alter  the  gust  response  of  an  air- 
plane from  that  obtained  by  considering  the  single  degree  of 
freedom  of  vertical  motion  only.  Results  can  be  expressed  in  a 
form  which  is  essentially  independent  of  the  Integral  scale  of 
turbulence.  It  appears  feasible,  therefore,  to  eliminate  the 
turbulence  scale  - which  has  been  the  subject  of  much  controversy  - 
as  one  of  the  parameters  in  the  design  for  gust  encounter. 

At  the  beginning  of  this  investigation  it  was  the  hope  that 
parametric  charts  could  be  derived  to  cover  a variety  of  aircraft 
configurations.  For  the  case  of  the  single  degree  of  freedom  of 
vertical  motion  only,  results  are  found  to  be  expressible  in  .terms 
of  a single  airplane  parameter,  the  mass  ratio  y . Addition  of 
the  degree  of  freedom  of  pitch  is  found  to  increase  the  number  of 
parameters  to  at  least  7-  This  large  number  essentially  precludes 
the  generation  of  parametric  charts.  On  the  other  hand,  the 
response  analysis  procedure  developed  herein  leads  to  fairly  simple 
response  equations,  as  given  in  Table  I.  Application  of  these 
equations  is  rather  simple  with  modern  computing  machines  and,  thus, 
results  for  various  configurations  can  be  established  almost  as 
quickly  as  using  parametric  charts. 

To  keep  results  tractable,  and  to  find  out  the  primary  signi- 
ficance of  including  pitch,  the  analysis  was  restricted  to  straight 
wings,  and  for  incompressible  flow.  The  analysis  is  easily  extended 
to  the  compressible  flow  ca3e.  The  only  change  would  be  tne  use  of 
different  definitions  of  cx  » si  » and  CQ  • The  equivalents  of 

these  quantities  for  compressible  flow  are  given  In  reference  2. 

The  extension  to  swept  wings,  or  to  delta  wing  planforms,  is  also 
quite  straightforward,  and  involves  mainly  the  consideration  of  the 
more  complex  geometry. 

The  analysis  herein  was  confined  to  two  line  loads  on  the  wing 
and  one  on  the  tail.  For  a more  detailed  consideration  of  a 
specific  airplane  configuration,  the  airloads  may  be  expressed  in 
terms  of  a more  detailed  grid  mash  than  used  herein.  For  example, 
to  account  for  spanwise  loading  effects  in  a more  precise  way  than 
used  in  this  study,  three  line  loads  on  each  half  span  at  each  of 
two  chord  positions  might  be  used  (for  symmetric  loading,  this 
choice  amounts  to  6 loads  as  compared  to  2 for  the  present  analysis). 
The  equations  that  would  result  wculd  be  Just  as  simple  as  those 
shown  in  Table  I,  and  simply  would  be  more  in  numoer  (one  for  each 
wing  load  assumed,  one  for  each  tail  load,  and  two  for  the  z and 
♦ motions).  It  Is  to  be  noted  that  the  type  of  analysis  used 
herein,  which  leads  to  the  rather  simple  equations  of  Table  I,  can 
also  be  used  to  include  flexibility  affectej  essentially,  one  more 
equation  would  be  added  for  each  flexible  mode  that  is  considered. 
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Figure  1,-  Influence  of  Various  Degrees  of  Freedom 


on  Response  Spectrum 
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(b)  Convolution  for  downwash 
ahead  of  line  load 


(c)  Convolution  for  downwash 
aft  of  line  load 


Figure  3.~  Geometry  for  Determination  of  Downwash  Due  to  Line  Load 


Figure  5--  Airplane  Transfer  Function  Check  for  One  Degree  of  Freedcr. 
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Pigujrs  6.-  Corapar 


Figure  7.-  Distribution  of  Response  Power  for  One  and  Two  Degrees  of  Freedom 
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Figure  8.-  K Curves  for  Case  III 


Figure  9.-  k Curves  fcr  Case  III 
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Figure  10. - K Curves  for  Cases  I,  II,  III  and  IV 
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Curves  for  Cases  with  Parameter  Variations 


Figure  13.-  k Curves  for  Cases  with  Parameter  Variations 


Effect  of  Tail  Downwash  on  Airplane  Transfer  Function 
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Comparison  of  kc  Curves  for  Airplanes  with  One  and  Two  Degrees  of  Freedom 


